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nonlinear  critical  stress  intensity  factor  (as  defined  here)  was  independent 
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of  the  size  of  the  specimen  and  was  constant  during  the  post-critical  crack 
growth. 


2.  Research  Objective 


To  develop  rational  methods  to  predict  fracture  toughness  of  fiber 
reinforced  concrete  is  the  primary  goal  of  this  investigation.  It  is  known 
that  the  addition  of  randomly  distributed,  short,  steel  fibers  substantially 
enhances  the  crack  propagation  resistance  of  concrete.  How  to  experimentally 
evaluate  this  critical  property  of  fiber  reinforced  concrete  and  how  to 
rationally  predict  it  are  the  objectives  of  this  research. 


3.  Status  of  the  Research 


A  method  is  developed  to  predict  crack  propagation  in  unreinforced 
matrix.  This  method  was  based  on  experimental  results  of  notched-beam 
specimens . 

The  closed-loop  notched-beam  tests  are  now  being  conducted  for  fiber 
reinforced  concrete.  The  theoretical  model  developed  for  plain  concrete  will 
be  extended  for  fiber  reinforced  concrete.  The  theoretical  model  proposed  is 
for  Mode  I  crack  propagation.  Attempts  are  being  made  to  extend  it  to  mixed¬ 
mode  loading. 
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A  MEASURE  FOR  THE  FRACTURE  TOUGHNESS  OF  CEMENT  BASED  MATERIALS* 


by 


Y.  S.  Jenq*  and  S.  P.  Shah^ 

Department  of  Civil  Engineering 
Northwestern  University 
Evanston,  Illinois,  60201 

Introduction 

It  is  frequently  reported  that  the  higher  the  strength  of  cement  based 
materials,  the  more  brittle  is  their  behavior.  It  could  be  useful  to 
quantitatively  express  the  degree  of  brittleness.  Many  attempts  [1-13]  have 
been  made  to  use  linear  elastic  fracture  mechanis  (LEFM)  to  quantitatively 
express  the  degree  of  brittleness.  For  example,  by  testing  notched  beams 
one  can  calculate  using  the  formulas  developed  from  LEFM,  a  quantity  called 
fracture  toughness  and  termed  *lc  from  the  measured  maximum  load  and  the 
initial  notch-length.  Unfortunately,  it  has  been  observed  that  thus 

calculated  is  dependent  on  the  dimension  of  the  beams.  Many  researchers 
have  attempted  to  analyze  this  size  dependency.  Such  approaches  are  usually 
quite  cumbersome  and  are  often  based  on  expensive  nonlinear  finite  element 
programs.  In  this  paper  a  direct  method  is  suggested  to  calculate  two  size- 
independent  fracture  toughness  parameters  from  the  experimental  results. 

The  method  was  developed  based  on  tests  on  notched-beams  of  different  mix 
proportions  and  different  sizes. 

C 

Values  of  KIc  <  critical  stress  intensity  factor)  and  CTODq  (critical 
crack  tip  opening  displacement)  are  suggested  which  appear  to  be  size 

C 

independent.  The  values  of  K?c  and  CT0Dc  (Fig.  1)  can  be  calculated  based 
on  the  peak  load,  initial  crack  length  (a0)  and  the  measured  crack  mouth 
opening  displacement  (CMOD).  The  calculation  is  based  on  the  observation 
that  nonlinear  crack  growth  [14-22]  occurs  and  that  the  elastic  displacement 

+For  presentation  at  the  Materials  Research  Society,  November  26-30,  1984. 

^Graduate  Research  Assistant,  Northwestern  University, 

Evanston,  Illinois  60201. 

^Professor  of  Civil  Engineering,  Northwestern  University, 

Evanston,  Illinois  60201. 


displacement  should  be  extracted  from  the  total  displacement  to  apply  linear 
elastic  fracture  mechanics. 

c 

Based  on  the  concept  of  constant  Kjc  and  CT0Dc,  theoretical  prediction 
of  load-CMOD  relationships,  size  effect  on  modulus  of  rupture  and 
conventional  critical  stress  intensity  factor  (Kjc)  were  compared  with  the 
experimental  results  and  good  correlations  were  observed. 

Inelastic  Displacement  and  Slow  Crack  Growth 

A  plot  of  load  vs.  crack  mouth  opening  (CMOD)  of  a  concrete  beam 
obtained  from  the  experiments  (Fig.  2)  can  be  used  to  describe  the  observed 
inelastic  displacement  and  slow  crack  growth  phenomena.  A  significant 
amount  of  inelastic  displacement  can  be  observd  when  the  specimens  are 
unloaded  immediately  after  the  peak  load.  The  total  crack  mouth  opening 
displacement  (CMOD^)  is  composed  of  a  sum  of  the  elastic  crack  mouth  opening 
displacement  without  slow  crack  growth  (CM0De),  inelastic  crack  mouth 

it  ^ 

opening  displacement  (CMOD  ),  and  the  elastic  crack  mouth  opening 
displacement  due  to  slow  crack  growth  (CMOD®)  (as  shown  in  Fig.  2).  The 
nonlinear  displacement  can  be  due  to  creep,  microcracking  and  slow  crack 
growth.  For  the  short-term  loading  considered  here,  it  is  assumed  that 
creep  effect  is  negligible.  The  phenomenon  of  slow  crack  growth  prior  to 
the  peak  load  has  long  been  noticed.  In  addition  to  the  geometrical 
interlock  and  microcracking  resulting  from  the  heterogeneous  nature  of 
concrete,  the  differences  in  the  stress  state  along  the  crack  front  can  also 
account  for  this  observed  slow  crack  growth.  It  has  been  noticed  that  the 
length  of  the  crack  front  at  the  surface  of  the  beams  is  longer  than  that  at 
the  interior.  Thus,  a  channel  shaped  slow  crack  growth  profile  was  observed 
by  Swartz  and  Go  [19].  Therefore,  it  is  difficult  to  precisely  measure  the 
length  of  the  crack  at  the  peak  load  using  optical  measurement. 

C 

Determination  of  Kjc  and  CT0Dc  Using  Three-Point  Bend  Tests 

It  is  clear  that  in  order  to  apply  LEFM,  the  elastic  crack  mouth 
displacement  (CM0De  =  CMOD^  -  CMOD  )  should  be  extracted  from  the  total 
displacement  (CMOD^).  Also,  to  overcome  the  difficulties  in  measuring  the 
exact  length  of  the  crack,  an  effective  crack  length  is  defined  in  this 
paper.  The  effective  crack  length  (a)  is  the  sum  of  the  initial  notch 
length  (ag)  plus  an  effective  crack  extension  (*e).  For  a  three-point 
bend  test  of  notched-beam  with  span-depth  ratio  of  four,  the  elastic  crack 


mouth  opening  displacement  can  be  expressed  by  an  empirical  formula  with 
accuracy  of  1%  error  [23]  as: 
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CM00e  =  6/-S--  V,  (§) 
b2  t  E'  1  b 


(1) 


where 


V,(4)  =  0.76  -  2.28A  +  3.87A2  -2.04A3  +  0,66 — / 
1  o  (1  -  A)2 

E'  -  Young’s  modulus  of  elasticity,  A  =  (■g-) 


P,  s,  a,  b,  t  are  indicated  in  Fig.  3. 


1 

I 
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For  a  given  measured  peak  load  (Pmax)*  initial  notch  lngth  (ap)  and  the 
measured  elastic  CM0De,  the  effective  crack  length  (a)  is  determined  so 
that  the  calculated  CM0De  is  equal  to  the  measured  CM0De.  Once  the 
effective  crack  length  is  determined,  then  one  can  calculate  using  LEFM  the 
values  of  and  CT0Dc.  The  k\c  can  be  calculated  as  [23] 


1 *5  Pmax  s 


F1 


in  which 


1  1.99  -  A(1  -  A)  (2.15  -  3.93A  +  2.7A2) 

/?  (1  +  2A)  (1  -  A) 2 


(2) 


The  value  of  CT00c  is  the  crack  opening  displacement  at  original  notch 
tip  and  can  be  calculated  using  LEFM  [26], 

Note  that  Kjc  is  not  the  same  as  the  conventional  critical  stress 
intensity  factor  CkIc)  because  calculation  of  Kjc  includes  effects  of  slow 
crack  growth  and  nonlinear  displacement.  To  calculate  Kjc  from  test 
results,  both  the  load  and  crack  mouth  opening  displacement  have  to  be 
monitored.  Unlike  the  conventional  values  which  were  found  to  be  size- 

dependent  [1-8,  10-13],  the  values  of  K3C  as  well  as  CT0Dc  were  found  to  be 
size-independent.  In  addition,  using  the  constant  values  of 
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Kjc  and  CTODc  it  was  possible  to  accurately  predict  the  load  vs.  CMOD  ^as 
well  as  load  vs.  load-point  deflection)  relationship  for  beams  tested  in 
this  investigation  as  well  as  by  other  researchers  [24,  25].  The 
theoretical  justification  for  the  validity  of  Kjc  is  given  in  Ref.  24. 

Test  Program  and  Experimental  Details 

Three  point  bend  test  was  used  to  verify  the  validity  of  the  proposed 
fracture  criteria.  Three  different  sizes  of  beams  were  used.  These  beams 
were  designated  as  large  (L),  medium  (M),  and  small  (S).  The  span  (s), 
depth  (b),  and  thickness  (t)  for  these  beams  were  respectively;  36  in. 
(914mm),  9  in.  (229mm),  and  3.375  in.  (85.7mm);  24  in.  (609mm),  6  in. 
(152mm),  2.25  in.  (57.2mm);  and  12  in.  (305mm),  3  in.  (76mm),  and  1.125  in. 
(28.6mm)  (Fig.  3).  Four  series  (Cl,  Ml,  M2,  PI)  of  different  mixes  were 
prepared  for  studying  the  effect  of  the  maximum  aggregate  size  and  water- 

.  C 

cement  ratio  on  the  fracture  toughness  Kjc  and  CT0Dc.  The  notches  were 
precast  for  Cl,  Ml,  M2  and  PI  series  and  were  saw  cut  for  M3  series.  All 
the  initial  notch-depth  was  equal  to  one-third  of  their  depth.  The  beams 
were  tested  after  curing  for  approximate  90  days  (Cl,  Ml,  M2,  and  PI)  and 
75  days  for  M3  series  in  an  environment  with  96£  relative  humidity  and  80°F 
(26.7°C).  The  mix  proportions  are  listed  in  Table  1. 

The  large  and  medium  size  specimens  were  tested  in  a  closed-loop  servo 
controlled  testing  machine  with  a  capacity  of  120  Kips  (534  KN)  while  the 
small  beams  were  tested  in  a  servo  controlled  closed-loop  testing  machine 
with  a  capacity  of  20  Kips  (89  KN).  All  notched  beams  were  tested  so  as  to 
maintain  a  constant  rate  of  increase  of  CMOD,  which  was  measured  by  a  clip 
gage.  The  unnotched  beams  were  tested  by  a  constant  rate  of  increase  of 
deflection.  The  peak  load  was  reached  in  about  10  minutes.  During  testing 
the  deflection  of  the  beam  was  measured  by  an  LVDT  located  at  about  1.5  to 
2  in.  (38  to  51mm)  away  from  the  notch  and  right  under  the  center  line  of 
the  beam.  The  surface  crack  length  was  monitored  by  a  microscope  with  an 
accuracy  of  0.001  in.  (0.0025mm).  The  overall  set-up  is  shown  in  Fig.  4. 
There  were  at  least  two  beams  for  each  size  and  mix  proportion;  a  total  of 
46  beams  were  tested  and  analyzed.  For  each  size  and  each  mix  proportion  at 
least  one  beam  was  loaded  cyclically. 

Discussion  of  Test  Results 


Some  of  the  test  results  are  summarized  in  Table  2.  Table  2  shows 
values  of  fracture  toughness  Kjc  and  CT0Dc  calculated  using  the  procedure 


c 

described  earlier.  It  can  be  seen  that  the  values  of  Kjc  as  suggested  here 
are  independent  of  the  size  of  the  beam.  This  was  not  true  for  the  values 
of  the  conventionally  calculated  using  the  initial  notch  length  ag  and 
the  maximum  load.  The  critical  crack  tip  opening  displacement  (CT0Dc)  were 
also  found  to  be  size-independent  (Table  2)  for  all  the  testing  series 
except  for  the  large  beams  of  Cl  series. 

c 

From  Table  2,  it  can  be  observed  (by  comparing  the  Klc  values  of  Cl 
series  vs.  Ml  series  or  M2  series  vs.  PI  series)  that  with  the  same  water- 
cement  ratio,  the  larger  the  size  of  the  aggregate  is  the  higher  are  the 

s  s 

values  of  fracture  toughness  Kjc.  Also,  by  comparing  Kjc  values 
associated  with  different  compressive  strength  (or  w/c  ratio)  of  the  same 
aggregate  size,  (Ml,  M2,  and  M3  series),  it  was  observed  that  values  of 

C 

Kjc  increase  as  the  values  of  compressive  strength  increase.  However,  the 
increase  of  Kjc  values  are  not  directly  proportional  to  the  increase  of 
compressive  strength  or  water/cement  ratio.  It  was  noticed  that  values  of 
CT0Dc  are  about  the  same  for  the  same  aggregate  size  with  different 
water/cement  ratio.  CT0Dc  increases  with  an  increase  of  aggregate  size. 

It  should  be  noted  that  unlike  classically  brittle  materials,  in  cement 
based  composites,  cracks  can  propagate  in  a  "stable"  manner  even  beyond  the 
peak  load.  This  post-critical  crack  propagation  (termed  strain  softening) 

C 

can  be  predicted  with  a  constant  Kjc  as  shown  in  Fig.  5.  The  details  are 
given  in  Refs.  24  and  25. 


Some  Predictions  Based  on  the  Proposed  Fracture  Criteria 

C  ... 

Using  the  proposed  fracture  criteria  Kjc  and  CT0Dc,  it  is  possible  to 

calculate  the  peak  load  for  any  given  geometry  [26].  To  demonstrate  this,  a 

theoretical  analysis  on  the  three-point  bend  test  was  conducted.  For  this 

analysis  it  was  assumed  that  concrete  and  cement  paste  had  the  same  Young's 

modulus  (3  x  10^  psi )  and  the  same  Kjc  value  (1200  psi  /Tn).  Partly 

based  on  the  test  results  reported  here  it  was  assumed  that  the  elastic 

critical  crack  tip  opening  displacements  (CT0Dc)  for  concrete  and  cement 

paste  were  respectively  0.0006  in.  and  0.0001  in.  From  the  proposed  model, 

the  peak  load,  and  hence  the  values  of  the  conventional  K.  can  be 

—  S  c 

calculated.  Fig.  6  shows  the  ratio  of  Kj  /<Ic  for  different  specimen  sizes 
with  notch-depth  ratios  ranging  from  0.1  to  0.8  and  depth-span  ratio  of 
four.  It  can  be  seen  that  "Kjc  reaches  its  maximum  values  at 
aQ/b  =  0.2  ~  0.3.  Such  a  trend  has  been  reported  from  the  test  data  by 
Carpinteri  [27]  (Fig.  7)  and  several  other  investigators  [1,  2,  8,  10]. 
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Fig.  8  gives  the  predicted  results  of  notch-sensitivity,  (i.e.,  the  ratio  of 
flexure  strength  of  notched-specimen  and  failure  flexure  strength  of 
unnotched  specimens),  for  concrete  and  cement  paste.  It  can  be  noticed  that 
cement  paste  is  more  notch-sensitive  than  concrete.  Similar  notch- 
sensitivity  results  were  also  reported  by  Shah  and  McGarry  [5]  (Fig.  9)  and 
other  investigators  [10,  27].  Size  effect  on  modulus  of  rupture  (MOR)  of 
concrete  is  shown  in  Fig.  10.  Note  that  MOR  decreased  with  the  increasing 
specimen  depth  as  has  been  commonly  observed  from  test  results  (Fig.  11). 

Conclusion 


1.  Nonlinear  crack  growth  effect  in  the  process  zone  and  inelastic 
behavior  of  concrete  need  to  be  considered  in  order  to  apply  LEFM  to 
concrete. 

2.  The  fracture  parameters  Kjc  (critical  stress  intensity  factor)  and 
CT0Dc  (critical  crack  tip  opening  displacement)  calculated  as  proposed 
here,  are  found  to  be  size-independent. 

C 

3.  Based  on  concept  of  constant  Kjc  and  CT0Dc,  the  reported  size  effct  on 

conventional  modulus  of  rupture  and  notch-sensitivity  can  be 

accurately  predicted. 

Acknowledgment 


The  research  reported  here  was  supported  by  the  Air  Force  Office  of 
Scientific  Research  [Lt.  Col.  Lawrence  D.  Hokanson,  Program  Manager] 
Grant  No.  AFOSR  -  82-0243. 


References 


1.  Kaplan,  M.  F.,  "Crack  Propagation  and  the  Fracture  of  Concrete," 
Journal  of  the  American  Concrete  Institute,  Vol.  58,  No.  5,  pp.  591- 
610,  November  1961. 

2.  Lott,  J.  1.,  and  C.  E.  Kesler,  "Crack  Propagation  in  Plain  Concrete," 
T.  &  A.  M.  Report,  No.  648,  University  of  Illinois  at  Urbana-Champaign 
August  1964. 

3.  Moavenzadeh,  F.,  and  R.  Kuquel,  "Fracture  of  Concrete,"  Journal  of 
Material,  JMLSA.Vol.  4,  No.  3,  September  1969,  pp.  497-519. 

4.  Naus,  0.  J.,  and  J.  L.  Lott,  "Fracture  Toughness  of  Portland  Cement 
Concretes,"  ACI  Journal,  June  1969,  pp.  481-489. 

5.  Shah,  S.  P.,  F.  J.  McGarry,  "Griffith  Fracture  Criterion  and  Concrete, 
Journal  of  Engineering  Division  (ASCE),  Dec.  1971,  pp.  1663-1675. 

6.  Walsh,  P.  F.,  "Fracture  of  Plain  Concrte,"  Indian  Concrete  Journal, 
Vol.  46,  No.  11,  469,  470,  476,  Nov.  1972. 

7.  Brown,  J.  H.,  "Measuring  the  Fracture  Toughness  of  Cement  Paste  and 
Mortar,"  Magazine  of  Concrete  Research,  Vol.  24,  No.  81,  1972,  pp.  185 
196. 

8.  Higgins,  D.  D.,  and  J.  E.  Bailey,  "Fracture  Measurements  on  Cement 
Paste,"  Journal  of  Material  Science,  Vol.  11,  1976,  pp.  1995-2003. 

9.  Mindess.  S.,  and  J.  S.  Nadeau,  "Effect  of  Notch  Width  on  K.  for 
mortar  and  Concrete,"  Cement  and  Concret  Research,  Vol.  6,  1C 

1976,  pp.  529-534. 

10.  Gjtfrv,  0.  E.,  S.  I.  Sorensen,  and  A.  Arnesen,  "Notch  Sensitivity  and 
Fracture  Toughness  of  Concrete,"  Cement  and  Concrete  Research,  Vol.  7, 

1977,  pp.  333-344. 

11.  Swartz,  S.  E.,  K.  K.  Hu,  and  G.  L.  Jones,  "Compliance  Monitoring  of 
Crack  Growth  in  Concrete,"  Journal  of  the  Eng.  Mech.  ASCE,  Vol.  104, 
No.  EM4,  August  1978,  pp.  789-800. 

12.  Swartz,  S.  E.,  K.  K.  Hu,  and  C.  M.  J.  Huang,  "Stress  Intensity  Factor 
for  Plain  Concrete  in  Bending  -  -  Prenotched  versus  Precracked  Beams," 
Experiment!  Mechanics,  Nov.  1982,  pp.  412-417. 

13.  Strange,  P.  C.,  and  A.  H.  Bryant,  "Experimental  Tests  on  Concrete 
Fracture,"  Technical  Notes,  Journal  of  Eng.  Mech.  Div.,  ASCE,  Vol.  105 
No.  EM2,  April  1979,  pp.  337-342. 

14.  Wecharatana,  M.,  and  S.  P.  Shah,  "Double  Torsion  Tests  for  Studying 
Slow  Crack  Growth  of  Portland  Cement  Mortar,"  Cement  and  Concrete 
Research,  Vol.  10,  1980,  pp.  833-844. 

15.  Wecharatana,  M.,  and  S.  P.  Shah,  "Slow  Crack  Growth  in  Cement 
Composites,"  Journal  of  Structural  Div.,  ASCE,  June  1982,  pp.  1100- 
1113. 

16.  Wecharatana,  M.,  and  S.  P.  Shah,  "Prediction  of  Nonlinear  Fracture 
Process  Zone  in  Concrete,"  Journal  of  EMD,  ASCE,  Vol.  109,  No.  5, 
October  1983,  pp.  1231-1246. 


12 


17.  Wecharatana,  M.,  and  S.  P.  Shah,  "A  Model  for  Predicting  Fracture 
Resistance  of  Fiber  Reinforced  Concrete,"  Cement  and  Concrete  Research, 
Vol.  13,  1983,  pp.  819-829. 

18.  Ballarini,  R.,  S.  P.  Shah,  and  L.  M.  Keer,  "Crack  Growth  in  Cement 
Based  Composites,"  Engineering  Fracture  Mechanis,  (To  appear). 

19.  Go,  C.  G.,  and  S.  E.  Swartz,  "Fracture  Toughness  Techniques  to  Predict 
Crack  Growth  and  Tensile  Failure  in  Concrete,"  Report  150,  College  of 
Engineering,  Kansas  State  University,  July  1983. 

20.  Hillerborg,  A.,  M.  Modeer,  and  P.  E.  Petersson,  "Analysis  of  Crack 
Formation  and  Crack  Growth  in  Concrete  by  Means  of  Fracture  Mechanics 
and  Finite  Elements,"  Cement  and  Concrete  Research,  Vol.  6,  1976,  pp. 
773-782. 

21.  Hillemier,  B.,  and  H.  M.  Hilsdorf,  "Fracture  Mechanics  on  Concrete 
Composites,"  Cement  and  Concrete  Research,  Vol.  7,  1971,  pp.  523-536. 

22.  Gopalaratnam,  V.,  S.,  and  S.  P.  Shah,  "Softening  Response  of  Plain 
Concrete  in  Direct  Tension,"  Accepted  for  publication,  ACI  Journal, 
(MS5652) . 

23.  Tada,  H.,  P.  C.  Paris,  and  G.  R.  Irwin,  The  Stress  Analysis  of  Cracks 
Handbook,  Del  Research  Corporation,  Hellertown,  Pennsylvania. 

24.  Jenq,  Y.  S.,  and  S.  P.  Shah,  "A  Fracture  Toughness  Criterion  For 
Concrete",  accepted  for  publication.  Engineering  Fracture  Mechanics. 


25.  Jenq,  Y.  S.,  and  S.  P.  Shah,  "Nonlinear  Fracture  Parameters  for  Cement 
Based  Composites:  Theory  and  Experiments,"  Proceedings  of  NATO  Workshop 
on  “The  Application  of  Fracture  Mechanics  to  Cementitious  Composites," 
edited  by  S.  P.  Shah,  pp.  213-253,  1984. 

26.  Jenq,  Y.  S.,  and  S.  P.  Shah,  (in  preparation) . 


27.  Carpi nteri.  A.,  "Static  and  Energetic  Fracture  Parameters  for  Rock  and 
Concrete,"  Materiaux  et  Constructions,  Vol.  14,  No.  81,  1981. 

28.  Wright,  P.  J.  A.,  Magazine  of  Concrete  Research,  Vol.  4.,  No.  11,  1952, 
pp.  67-76. 


Load,  P  (KN) 


1200  pal  fin  (1. 319MN 


Wright  [28] 


END 

FILMED 

2-85 

DTIC 


